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ABSTRACT 

Simultaneous spectral differential imaging is a high contrast technique by 
which subtraction of simultaneous images reduces noise from atmospheric speck- 
les and optical aberrations. Small non-common wave front errors between chan- 
nels can seriously degrade its performance. We present a new concept, a multi- 
color detector assembly (MCDA), which can eliminate this problem. The device 
consists of an infrared detector and a microlens array onto the flat side of which a 
checkerboard pattern of narrow-band micro-filters is deposited, each micro-filter 
coinciding with a microlens. Practical considerations for successful implemen- 
tation of the technique are mentioned. Numerical simulations predict a noise 
attenuation of 10~ 3 at 0.5" for a 10 5 seconds integration on a m# = 5 star of 
Strehl ratio 0.9 taken with an 8-m telescope. This reaches a contrast of 10 -7 at 
an angular distance of 0.5" from the center of the star image. 

Subject headings: instrumentation: adaptive optics - instrumentation: detectors 
- planetary systems - stars: imaging 



1. Introduction and Motivation 

Direct imaging and spectroscopy of exoplanets are crucial to constrain their masses, 
study their atmospheres and, ultimately, find telltale signatures of life. The task is dauntingly 
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difficult. The planet image is hugely fainter than the background from the brilliant stellar 
image. The most obvious problem is photon noise in the bright star point spread function 
(PSF). A remedy is coronagraphy (Lyot & Marshall 1933; Roddier & Roddier 1997; Rouan 
et al. 2000; Baudoz et al. 2000; Abe et al. 2001; Kasdin et al. 2003), which attenuates the 
coherent diffraction pattern of the on-axis PSF. At high Strehl ratios the diffraction pattern 
dominates much of the PSF and coronagraphs can achieve strong attenuation (Kilston et al. 
2002). With the relatively modest Strehl ratios achieved with current adaptive optics (AO) 
systems on large ground based telescopes, the efficiency of coronagraphy is limited by the 
presence of a rapidly changing halo of speckles generated by atmospheric phase distortion. 
This increases the local PSF variance by orders of magnitude above that caused by photon 
noise (Racine et al. 1999). In addition, aberrations in the optical train, unfiltered by the AO 
system, produce speckles of much longer lifetime, hence much brighter in a long-exposure 
image, than atmospheric speckles. This noise can be dominant for integrations longer than 
a few seconds (Marois et al. 2003a, 2004). Quasi-static speckles also limit the high-contrast 
performance of the Hubble Space Telescope (Schneider & Silverstone 2003), in the form of 
the so-called "breathing" problem. All this underscores the importance of PSF calibration 
for high contrast imaging, which is the topic of this paper. 

Simultaneous spectral differential imaging (Smith 1987; Racine et al. 1999; Marois et 
al. 2000a; Sparks & Ford 2002; Biller et al. 2004) is such a calibration technique. Images 
are acquired simultaneously at adjacent wavelengths in a spectral range where the planetary 
and stellar spectra differ appreciably. Judicious image subtraction removes the stellar image 
and conserves that of any companion. A useful spectral feature for this purpose is the sharp 
methane absorption band head at 1.6 /im (Rosenthal et al. 1996) found in relatively cold 
atmospheres (Burgasser et al. 2002) such as those of cold T type brown dwarfs and Jovian 
planets. 

The TRIDENT camera (Marois et al. 2000b, 2003a,b,c, 2004), a 3-channel (1.580, 1.625 
and 1.680 /im, 1% bandwidth) differential imager, has been extensively used with the AO 
system of the 3.6-m Canada- France-Hawaii Telescope (Rigaut et al. 1998). From the data, 
one can determine the level of noise attenuation achieved. The noise attenuation factor 
AN/N is defined as the median over an annulus of the pixel- by-pixel ratio of the absolute 
value of the intensity in the difference image to its absolute value in the original image 
after subtraction of an azimuthally averaged profile. The noise structure of the PSF at one 
wavelength is dominated by static aberrations and is well approximated with ~ 130 nm of 
static phase screen following a power-law distribution P oc u a with a = —2.7. This level of 
aberration is consistent with that measured (110 nm RMS) for the static component of the 
PUEO AO system (Rigaut et al. 1998). Fig. 1 illustrates typical noise attenuations obtained 
with TRIDENT. At subarcsecond separations the noise residuals in a difference image are 
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~ 10 times larger than the sum of the photon and read noises. Results from numerical 
simulations indicate that ~ A/20 of non-common path aberrations can explain TRIDENT 
limited performance. A factor of ~ 4 improvement is obtained by subtraction of the residuals 
of a reference star observed shortly after and under the same conditions as the science target. 
This indicates that differential aberrations between channels contribute a large fraction of 
the residual noise in a difference of simultaneous images. It also appears that some of the 
noise after subtraction of the reference star image is caused by slight changes in optical 
aberrations between the target and the reference star observations. 

In this paper, in light of the non-common path limitations, a new differential imager 
concept is introduced that should improve detection limits by a large factor. Simulations 
are made to estimate the performances and the expected detection limits with 8-m class 
telescopes and high order AO systems. 

2. The Multicolor Detector Assembly Concept 

A solution to the non-common path problem is to perform the spectral separation in the 
image plane of a single optical path imager with a Multicolor Detector Assembly (MCDA). 
A three- wavelength MCDA is illustrated in Fig. 2. A micro lens array located at the focal 
plane is used to concentrate the light of each PSF spatial sample into well separated detector 
pixels (each microlens covers a 5 x 5 pixel cluster) in order to minimize interpixel crosstalk, 
as discussed below. A mask on both faces of the microlens array blocks the light which may 
scatter between microlenses. Each microlens in a 2 x 2 microlens cluster is covered with 
one of three different micro-filters, two microlenses having the same filter for signal-to-noise 
optimization. The fact that only a small central area of each micro filter is used relaxes 
the requirement for edge-to-edge filter uniformity. A detector read-out yields four critically 
sampled "checkerboard" images of the PSF, one at each wavelength, that can be interpolated 
to construct a filled image. 

The attenuation of noise structures allowed by this device, in the absence of photon, read 
and flat field noises, can theoretically be arbitrarily large if the noise structures - speckles - 
are critically sampled at each wavelength over a PSF image whose extent is unlimited. This 
is so because, by the sampling theorem for bandwidth limited images, it is then possible to 
perfectly interpolate the pixel intensity at a given wavelength for any pixel which is "blind" 
at that wavelength. Complete attenuation is of course not achievable in practice because the 
idealized conditions posed above cannot be realized. Flat field errors will probably be the 
dominant limitation at high flux levels. The MCDA itself may entail limiting factors such 
as the two mentioned below. 
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The microlens array is required to minimize signal crosstalk between PSF sampling 
points. With current infrared detectors, a pixel can leak as much as 40% of its signal to its 
adjacent pixels (Finger et al. 1998). Diffraction by the microlenses can also be a significant 
source of crosstalk. Crosstalk limits the PSF noise attenuation as follows. The distance 
of a speckle to the PSF center is proportional to wavelength. PSFs obtained at different 
wavelengths must be re-scaled before being subtracted. The leak of a speckle signal sampled 
at one wavelength to a sampling point for another wavelength no longer coincides spatially 
with itself after re-scaling. Subtraction of the re-scaled images leaves a residual whose 
intensity is equal to the crosstalk coefficient X times the difference in the speckle intensity 
between the radially shifted samples. The average residual noise after image subtraction 
then increases from zero at the PSF center, where the wavelength shift is nil, to X times its 
initial level when the radial shift between wavelengths exceeds the FWHM (~ A/D) of the 
speckles, i.e. at an angular offset 9 > A 2 /(AA x D). At larger offsets, the crosstalk induced 
noise saturates at an amplitude of X times the speckle intensity and the limiting attenuation 
factor becomes equal to the crosstalk level X. Numerical simulations (Fig. 3) confirm this 
analysis. It is expected that sample separations ~ 5x 18.5 /iin pixels (to reduce charge leaks 
between samples to < 0.001) and microlens f-ratios < 17 (to limit the spillover of a pupil 
diffraction pattern onto an adjacent pupil image) will allow X < 0.001 at A ~ 1.6 /im. 1 To 
minimize the sample to sample crosstalk and to block gaps between microlenses, two masks 
are introduced above and under the microlens array. 

PSF decorrelation between wavelengths can also occur with an MCDA if the optical 
path is chromatic. Refractive optics produce wave fronts of different dimensions at different 
wavelengths because of the wavelength dependence of the index of refraction. Atmospheric 
dispersion also decorrelates the PSFs noise structure because the wave fronts for different 
wavelengths follow different optical paths through the atmosphere and instrument. This 
leads to non-common path errors since the footprints on optical surfaces of wave fronts at 
different wavelengths will be sheared. Simulations indicate that a ~1/1000 wave front shear 
limits attenuation to ~10~ 3 . The effects of chromatism and atmospheric refraction can be 
mitigated by using reflective optics and atmospheric dispersion compensators. 



1 For a circular aperture, the mean intensity diffracted at radius r is ~ 0.04 [r/ (A x //D)] J for unit 
central intensity. For diffraction crosstalk to be < 0.001 at a wavelength of 1.6 /im and at a distance of 
5x18.5 fxm requires //D < 17. 
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3. Simulated Performance with an MCDA 

Numerical simulations were done to estimate the PSF noise attenuation achievable with 
a three-wavelength (Ai = 1.52 /zni, A2 = 1.58 fim and A3 = 1.64 fim) MCDA in quasi-realistic 
conditions. Simulations are for a 10 5 s total integration divided in 60 s exposures on an 
m H = 5 star observed with an 8-m telescope and 4k actuator AO system of 20% transmission 
and 2% bandpass. AO phase filtering is simulated by choosing a flat power spectrum for 
spatial frequencies lower than twice the interactuator spacing. Static aberrations are included 
using a power- law power spectrum with a = —2.7 to obtain 100 nm RMS total aberrations, 
of which 15 nm RMS are not corrected by the AO system due to non-common path optics. 
The three images are the average of 500 atmospheric speckle realizations for 0.5" seeing at 
0.5 /im. AO filtering of both the static and Kolmogorov atmospheric phase screens yields 
PSFs with Strehl ratios of 0.96 at 1.60 /xm. These also include a typical level of read noise 
of 5 e- per pixel per 60 s exposure, photon noise, random flat field errors of 10~ 4 , sky noise 
and sample to sample crosstalk of 10~ 3 . PSF evolution with wavelength is included and 
partly compensated by appropriate combination of the multi-wavelength images (Marois et 
al. 2000a; Marois 2004). For the sake of completeness, simulations were done both with 
and without an optimized Lyot coronagraph. For the coronagraph simulation, the 15 nm 
RMS unfiltered static aberration is split in two uncorrelated 10 nm RMS aberrations that 
are respectively included before and after the coronagraph. The results are shown in Fig. 4. 
These lead to the detection limits (Fig. 5), checked with off-axis companion images, which are 
compared to the contrast required for the detection of 2.5 Mj up exoplanets around 0.5 Gyr 
M0 stars (Baraffe et al. 2003). When combined with a MCDA, the coronagraph does ~ 1 
magnitude better. The coronagraph has effectively attenuated the PSF structure by an order 
of magnitude inside a 0.5" radius, thus reducing the flat field accuracy needed to achieve 
AH = 17.5. 



4. Summary and Conclusion 

Experience with the TRIDENT camera has shown that small non-common path aber- 
rations in multiple optical channel instruments seriously limit PSF noise attenuation hence 
faint companion detection. The MCDA concept described here can overcome this problem. 
Numerical simulations predict that, with an MCDA, 2.5 Mj up exoplanets around 0.5 Gyr 
M0 stars should be detectable in 10 5 seconds with an 8-m telescope and a 4k actuator AO 
system. 

Planet finders based on MCDAs would enable detection and characterization of exo- 
planets with present large telescopes and high order AO systems. The technique is relatively 
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simple to implement and would complement coronagraph projects, or telescopes like HST, 
where it could calibrate evolving noise structure. An MCDA combined with a coronagraph 
would be a powerful tool for detecting ozone in the mid-IR (Des Marais et al. 2002), water 
in the near-IR, molecular oxygen in the red and, eventually, the chlorophyll "red edge" with 
the Terrestrial Planet Finder. An optimized MCDA-type polarimeter could study reflected 
star light from exoplanets (Saar & Seager 2003), or disks and zodiacal dust. Instead of 
micro-filters, micro-polarizers would be used to produce simultaneous polarized images with 
a single optical channel. 

This work was supported in part through grants from NSERC, Canada and from 
FQRNT, Quebec. 
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Fig. 1. — Median of the pixel by pixel ratio of residual noise to initial noise as a function 
of the distance from the PSF center for (a) the difference of two simultaneous TRIDENT 
images (solid line), (b) the difference between differences of simultaneous TRIDENT images 
of two stars observed sequentially (dashed line), (c) a difference image whose noise would only 
consist of the speckle, read, photon and flat field noises of the TRIDENT images (dotted 
line) and (d) a simulated difference image with the same noises plus an RMS wave front 
difference between channels of A/20 (dot-dash line). 
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Fig. 2. — Concept of a three wavelength differential imager with an MCDA showing the 
whole assembly (upper right) and the micro-filter array (upper left). Here, each microlens 
covers a 5 x 5 pixel cluster. A PSF and one of its checker board monochromatic samples are 
shown at the bottom. See text for details. 
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Fig. 3. — Effect of sample to sample crosstalk in an MCDA from numerical simulations. This 
shows the noise attenuation achieved by the subtraction of two PSF samples as a function 
of speckle separation with wavelength, expressed in A/D units, when attenuation is only 
limited by 0.001 (dotted line), 0.01 (dashed line) and 0.1 (solid line) of the sample to sample 
crosstalk. 
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Fig. 4. — PSF noise attenuation around an m H = 5 star observed for 10 5 second with an 
8-m telescope, a high order AO system delivering a Strehl ratio of 0.91, a three-wavelength 
MCDA and 20% global optical transmission without (solid line) and with (dashed line) a 
Lyot coronagraph. The latter features a gaussian mask with a FWHM of 15 X/D and a 60% 
transmissive Lyot mask. 
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Fig. 5. — 6(7 detection limit with (three dots dashed line) and without (dot dashed line) 
a Lyot coronagraph after multicolor speckle suppression. The dashed ans solid lines show 
respectively the noise of a single PSF with and without a Lyot coronagraph. The dotted 
line shows the magnitude difference for a 0.5 Gyr / 2.5 Mj up (from Baraffe et al. (2003)) 
companion around a MO star at 10 pc, including the effect of stellar flux reflection by the 
companion at small separations. 



